Interface Depinning in the Absence of External Driving Force 
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We study the pinning-depinning phase transition of interfaces in the quenched Kardar-Parisi- 
Zhang model as the external driving force F goes towards zero. For a fixed value of the driving force 
we induce depinning by increasing the nonlinear term coefficient A, which is related to lateral growth, 
up to a critical threshold. We focus on the case in which there is no external force applied {F = 0) 
and find that, contrary to a simple scaling prediction, there is a finite value of A that makes the 
interface to become depinned. The critical exponents at the transition are consistent with directed 
percolation depinning. Our results are relevant for paper wetting experiments, in which an interface 
gets moving with no external driving force. 
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The dynamics of random interfaces in the presence of 
noise is an interesting example of critical phenomena and 
generic scale-free behaviour in systems far from equilib- 
rium. In the case of surface growth dominated by ther- 
mal fluctuations, the Kardar-Parisi-Zhang (KPZ) equa- 
tion has been very much studied for it represents a 
whole universality class of growth, which includes many 
well-known discrete computer models . In many exper- 
imental situations, however, interface motion is affected 
by the existence of random pinning forces (see |^] and ref- 
erences therein). In this case, the simplest way to model 
interface roughening is to replace the noise term 77 (x, t) 
in KPZ by a quenched disorder 7y(x, /i). 



dh 
'dt 



vV^h + X(yhf + F + r;(x, h), 



(1) 



which is often referred to as the quenched Kardar-Parisi- 
Zhang (QKPZ) equation. The first term on the right- 
hand side describes the smoothening effect of surface 
tension, F is the driving force that pushes the interface 
through the disorder, and the term A(V/i)^ comes from 
lateral growth and represents the nonlinear most relevant 
correction. The quenched disorder has short-range corre- 
lations (77(x, /i)r;(x', h')) = S{x — x')A(/i — h'), where the 
correlator A(w) is a very rapidly decreasing function of |m| 
and is the term actually responsible for the pinning of the 
interface. This equation is expected to describe interface 
roughening in many disordered systems, including the 
non-equilibrium dynamics of magnetic domain walls in 
disordered materials |^-^ , an elastic chain in a quenched 
disorder , fracture cracks propagation ||] , etc) . Its ap- 
plicability to describing fluid-fluid displacement in porous 
media might be less justified though 

The QKPZ model described by (|l|) exhibits a continu- 
ous phase transition at a certain critical value, Fc, of the 
external driving force F. For F larger than Fc, the inter- 
face moves with a finite velocity. However, the interface 
remains pinned by the disorder for F < Fc- The critical 
point F = Fc is known as depinning transition. The in- 
terface velocity scales as f ~ (i^ — Fc)^ near and above 
the transition and plays the role of an order parameter. 

The value of the critical force depends on the param- 



eters of the model, in particular depends on the value 
of the coefficient A of the nonlinear term. Therefore, by 
keeping constant the rest of the equation parameters one 
can find a critical line Fc = /(A) separating the pinned 
from the depinned phase. Alternatively, we can see this 
critical line the other way around and let Ac = f~^{F) be 
the critical value of the KPZ nonlinearity above which the 
interface gets depinned. The driving force F favors the 
advance of the interface and thus, the lower the driving 
force is, the larger the critical value Ac of the nonlinear- 
ity that is needed in order to get the interface depinned. 
Indeed one would expect that as — > depinning be- 
comes more and more difficult until eventually, at F = 0, 
the threshold Ac — > 00 and depinning becomes impossi- 
ble. This intuitive picture can be justified by means of 
a simple scaling argument as follows. Consider a typical 
region of size I pinned by the disorder. Eq.([l]) applied to 
that region reads 

+ Xh^r^ + F- A(0)i/2;-d/2 ^ Q 

If one supposes that the nonlinear term dominates over 
the diffusion, the interface remains pinned whenever 
\a?l~'^ ^ A(0)-^/^Z~'^/^, where a is the lattice spacing 
in the growth direction. This defines a characteristic 
length, Ic = [A2aVA(0)]i/('^-'*), such that for I < Ic 
the interface gets pinned. Now to estimate the criti- 
cal force that is necessary to depin a region of typical 
size Ic, one equates the force term with the disorder 
in Eq.(||) to get to an expression for the critical line, 
Fc - A(0)2/(4-rf)(Aa2)-''/(4-d). Inverting the latter, one 
finds 

A(0)2/<i 



Ac 



_p-{A-d)/d 



io|. 



(3) 



In 



for the critical line of the depinning transition 
1 + 1 dimensions for instance, Eq.(|^) predicts a diverging 
Ac ^ as F -> 0. 

In this Letter, we show that, contrary to this scaling 
picture, there is always a finite critical value Ac of the 
KPZ nonlinearity such that the interface gets depinned 
even for _F = 0. Our conclusions are based upon numeri- 
cal integration of Eq.(|^) in d — 1. We numerically calcu- 
late the critical line and find that Xc{F = 0) = 3.60±0.01 
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(in natural units) for the QKPZ equation. Our results 
support the somehow counterintuitive conclusion that an 
interface can get depinned in absence of external driving 
force by the solely effect of nonlinearities. 

In order to numerically integrate Eq.(|l|), the equation 
parameters can easily be rescaled to have only two inde- 
pendent tuning parameters- namely, the nonlinear KPZ 
coefficient A and the driving force F. We have used 
a standard finite-differences scheme for integrating the 
QKPZ equation given (in natural units) by 

h{i,t + At) ^ h{i,t) + At F + At r][i,h{i,t)] (4) 
+ At [h{i + l,t) + h{i-l,t) -2h{i,t)] 
,_^ hii + l,t)-h{i-l,t) ^^ 

where the lattice spacing has been set to unity. We start 
our simulation from a flat initial condition h{x,0) — 
and periodic boundary conditions, i.e. /i(0,t) = h{L,t) 
and h{L + l,t) = h{l,t), are imposed on the inter- 
face. h{i,t) stands for the integer part of h{i,t), and the 
quenched disorder is Gaussian distributed and has corre- 
lations {r]{i,h)r]{j,h')) = SijSf^^,. Simulations with dif- 
ferent time steps were carried out, and the scheme proved 
to be stable and well behaved for a time step A = 0.01 (or 
smaller) for the range of tuning parameters simulated. 

We carried out simulations in systems of size L = 
128, 256, • ■ • , 8192. For each value of the of the nonlin- 
ear coefficient A we computed the critical value of force 
needed to get the interface depinned. Our results are 
summarized in Fig. 1. As expected, we find that as the 
driving force is smaller the critical value Ac of the non- 
linear coefficient required in order to depin the interface 
becomes larger. However, as anticipated above, the crit- 
ical point Ac always remains finite, even for F = 0. At 
a purely phenomenological level, we find that the critical 
line can be fitted very nicely by 
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where the constants 6i = 4.31 ± 0.04 and &2 = 0.81 ± 0.03 
(see Fig. 1). To our knowledge this is the first formula 
for the critical line and demands theoretical explanation. 

In the following we focus on the case in which no 
external driving, _F = 0, pushes the interface and de- 
pinning is due solely to nonlinear lateral growth. We 
have studied the critical behaviour in the vicinity of 
Ac(-F = 0) = 3.60 ± 0.01 in order to address the prob- 
lem of the nature of the critical point. Firstly, we have 
computed the scaling behaviour of the stationary inter- 
face velocity at = as the transition is approached. In 
Figure 2 (inset) we plot v vs. A for F = and a system 
of size L — 8192 showing that the transition is continu- 
ous. The critical behaviour of the order parameter v is 
shown in Figure 2. We find that close to the depinning 
threshold the interface velocity scales as u ~ (A — Ac)^ 
with a critical exponent 9 = 0.635 ± 0.007. 

The depinning mechanism for F ~ is the following. 
Starting from a flat initial condition h{x,t ^ 0) = all 



the terms in Eq.(|l|) are zero except for the disorder. At 
time t — the quenched random term ri[x,h) generates 
inhomogeneities in the front, which in turn produce a 
finite value of (V/i)^. For small values of A this inhomo- 
geneities smear out and the interface gets pinned by the 
disorder at one of the infinite pinning paths. However, for 
A > Ac these initial inhomogeneities are effectively am- 
plified by the nonlinearity and the interface gets moving 
with a finite velocity. 

As occurs in the standard case of depinning at a thresh- 
old value of the driving force F = Fc, the depinned phase 
is rough and belongs to the universality class of KPZ. 
This can be seen by studying the scaling behaviour of 
the the global width W{L, t) = [{h{x, tf) - {h{x, t))^]^/^^ 
where the average is over all x and different realizations 
of disorder . We obtain that the global width scales 
as 



W{L,t) 



tP 



if 
if 



(6) 



with a time exponent (3 = 0.33 ± 0.01 and a roughness 
exponent a — 0.50 ± 0.01 in agreement with the KPZ 
class of growth. 

However, when approaching the depinning transition 
from above, e = (A — Ac) /Ac 0+, the scaling of the 
global width is affected by the existence of a diverging 
correlation length ^ ^ e""^. This is the typical size of the 
fluctuations of the majority phase, i.e. the characteristic 
size of connected regions formed by pinned sites. As we 
show in Figure 3, the global width (and similarly does the 
local width) displays a crossover from ^ t^ '^ to KPZ-like 
behaviour ^ t'^-^^ . More precisely, one can see in Fig. 3 
that the width approximately behaves as 



W{t,e) 



if 
if 



t < tc 
t > tc 



(7) 



where Kc, in view of the dependence of the curves on e, 
must be very small. These two regimes are separated by 
a crossover time tc that depends on e. Indeed, following 
Kertesz and Wolf , near a roughening phase transition 
one expects the crossover time to scale with the distance 
to the threshold as tc ^ ^ where 7 = zi^. Direct 
examination of Figure 3 immediately suggests the scaling 
ansatz 



W(t,e) ~t^'=-e-«5(t/tc). 



(8) 



which is characteristic of systems close to a roughening 
transition ||l3-n5|. The scaling function is given by 
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y/3c-/3fcp. if M < 1 

const. if u 3> 1 



and the scaling relation 

Kc + K ^ iPc- Pkpzh 



(9) 



(10) 



among critical exponents must be fulfilled so that both 
regimes match. 

In Figure 3 (inset) we show a data collapse of 
t^^'"''e'^W{t,e) vs. e^t. A good data collapse is obtained 
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for the exponents Pkpz — 0.3, k — 0.57 and 7 — 1.57, the 
error in estimating these exponents being of about 10%. 
From the scahng relation ( p^ ) one also gets Pc = 0.73 in 
good agreement with our previous estimate. 

The value of the critical exponents is consistent with 
those of the DPD model (l^Jl^ just above the transition 
[p^P]. We thus conclude that the lateral growth driven 
depinning point at = and A = Ac also belongs to the 
universality class of DPD. 

Our results indicate that in the absence of any external 
driving field an interface can get depinned by increasing 
the nonlinear term A in up to its critical value. From 
the experimental point of view, this implies that, assum- 
ing the parameter A is tunable in the laboratory, an inter- 
face could become depinned even with when no external 
driving force is applied. In the following we discuss the 
role of anisotropy of the background random medium in 
generating the KPZ term A(V/i)^, and how this mecha- 
nism can be used to rise the value of A in experiments by 
increasing the degree of disorder anisotropy. 

The QKPZ equation for A = is known as the 
quenched Edwards- Wilkinson (QEW) equation and has 
been much studied in recent years. The critical expo- 
nents at the depinning transition have been well de- 
termined by several authors |18-2lj^. In 1+1 dimen- 
sions one finds a ^ 1.25 and f3 ~ 0.85 at the thresh- 
old F = Fc and a = 1/2 and (3 = 1/4 in the moving 
phase for F ^ F^., where the disorder ri{x,h) can be re- 
placed by r]{x, vt) and the exponents of the EW univer- 
sality class are recovered. The QEW equation arises 
naturally as the Langevin equation for the Hamiltonian 
H = J (ix[^l -I- (V/i)2 -I- V{x, h)] describing the elastic 
energy of an interface in a disordered potential l^(x, y) 
The term A(V/i)^ cannot be deduced as a variation 
of any Hamiltonian and is added as the most relevant 
nonlinear correction Geometrically, it accounts for 
growth in a direction locally normal to the interface and 
is referred to as nonlinear lateral growth term. 

In the past the physical origin of the KPZ nonlinear- 
ity in interface depinning has been found to be related 
to two distinct mechanisms for different models [Q. On 
the one hand, in the spirit of the original work of KPZ 
[0, the A term can have a purely kinematic origin, so 
that A oc w |l^j2^. In this case, the term A(V/i)^ goes 
to zero at the depinning transition, F — Fc, and the 
system thus belongs to the QEW universality class. On 
the other hand, there are models |2^ for which A re- 
mains finite at the transition |2^]. These models have 
exponents that correspond to the DPD universality class 



I JIS 24 1 . Tang, Kardar and Dhar |25| have shown that this 
finite A term can arise in some models because of an un- 
derlying anisotropy in the random medium, i.e. models 
that have a growth direction determined by the random 
medium. A further numerical step on this direction has 
recently been achieved by Park, Kim and Kim ||2^ by 
studying a model with an anisotropic disorder correla- 
tor. The effect of anisotropy on real experiments has 
also been successfully tested by Albert et. al. |2^. Ex- 
periments on fluid flow in a random medium formed by 



packed glass beads |2^ are now known to belong to the 
isotropic QEW universality class However, the scal- 
ing exponents obtained for paper wetting ||l^,^ are in 
excellent agreement with the prediction of the anisotropic 
DPD universality class. In paper wetting experiments a 
sheet of paper is vertically suspended over a reservoir of 
liquid (usually black ink) . The fluid then wets the paper 
and the interface between wet and dry phases rises until 
it eventually stops. The interface grows upwards because 
of capillarity forces in the paper pores. Notice that there 
is none external driving force. The anisotropic paper fi- 
bre distribution determines the local capillarity forces. 
Disorder in these systems is thus highly anisotropic. We 
believe that this system is an excellent example of depin- 
ning driven solely by the nonlinear lateral growth term. 

In summary, we have studied the QKPZ equation fo- 
cusing on the case in which there is no external driving 
force {F = 0). We have shown that there exists a de- 
pinning transition for a finite value of the KPZ coeffi- 
cient A = Ac(-F' = 0) and that transition belongs to the 
DPD universality class. A finite value of the nonlinear 
coefficient A appears in systems with anisotropic disor- 
der like for instance in paper wetting experiments. In 
this system there is no external driving force and depin- 
ning occurs due to local capillarity forces, which drive 
the interface through the anisotropic lateral growth term 
A(Vft.)^. We conclude that by varying the anisotropy 
degree of the corresponding random medium in other ex- 
perimental systems, depinning is possible even with no 
external driving. 
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FIG. 1. Critical line Ac = f{F) for the QKPZ equation. 
Symbols are points obtained from numerical simulations in a 
system of size L = 1024. The line is a fit according to Eq.(5). 
Note that Ac remains finite, even at _F = 0. 
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FIG. 2. Interface velocity vs. coefficient A for the QKPZ 
equation at _F = (inset) close to the threshold Ac(-F = 0). 
The critical behaviour of the velocity v ~ (A — Ac)'' is shown in 
the main panel. A straight line is found for Ac = 3.60 ± 0.01 
and the slope corresponds to the velocity critical exponent 
V = 0.635 ± 0.007. 




FIG. 3. In main panel we plot the global width for differ- 
ent distances (as shown) e = (A — Ac)/Ac to the threshold for 
_F = in a system of size L — 8192. The crossover from t"'^ to 

j.0.3 



t occurs at times that scale as tc 



with the distance 



to the threshold. Inset shows a data collapse according to 
Eq.(8) of the sets shown in the main panel. A good collapse 
is found for the exponents f3kpz ~ 0.3, k — 0.57 and 7 — 1.57. 
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